Multi-step forecasting of dissolved oxygen in River Ganga based on CEEMDAN-AdaBoost-BiLSTM-LSTM model

Accurate prediction of Dissolved Oxygen (DO) is an integral part of water resource management. This study proposes a novel approach combining Complete Ensemble Empirical Mode Decomposition with Adaptive Noise (CEEMDAN) with AdaBoost and deep learning for multi-step forecasting of DO. CEEMDAN generates Intrinsic Mode Functions (IMFs) with different frequencies, capturing non-linear and non-stationary characteristics of the data. The high-frequency and medium-frequency IMFs, characterized by complex patterns and frequent changes over time, are predicted using Adaboost with Bidirectional Long Short-Term Memory (BiLSTM) as the base estimator. The low-frequency IMFs, characterized by relatively simple patterns, are predicted using standalone Long Short-Term Memory (LSTM). The proposed CEEMDAN-AdaBoost-BiLSTM-LSTM model is tested on data from ten stations of river Ganga. We compare the results with six models without decomposition and four models utilizing decomposition. Experimental results show that using a tailored prediction technique based on each IMF’s distinctive features leads to more accurate forecasts. CEEMDAN-AdaBoost-BiLSTM-LSTM outperforms CEEMDAN-BiLSTM with an average improvement of 25.458% for RMSE and 37.390% for MAE. Compared with CEEMDAN-AdaBoost-BiLSTM, an average improvement of 20.779% for RMSE and 28.921% for MAE is observed. Diebold-Mariano test and t-test suggest a statistically significant difference in performance between the proposed and compared models.

The accuracy of a artificial intelligence models can be greatly enhanced by boosting and can aid in improving generalization and reducing overfitting.In the study by Aldrees et al. 13 , boosting and bagging ensemble models are compared with individual network-based and tree-based models for making predictions.El Bilali et al. 14 discussed the applicability of using AdaBoost and Random Forest for predicting groundwater quality.
Often water quality time series data exhibits large fluctuations, non-linearity, and non-stationary characteristics 15 which makes the task of forecasting challenging.To address this, pre-processing the data using signal decomposition techniques can be useful in transforming the complex dataset into relatively simple subseries.Sahoo et al. 16 , a uses Fourier Transform combined with LSTM to predict daily suspended sediment load.Huang et al. 17 introduced Empirical Mode Decomposition (EMD) which is a signal processing technique that analyzes non-stationary signals by breaking them into Intrinsic Mode Functions (IMFs) according to their local characteristics.Wu et al. 18 proposed the Ensemble Empirical Mode Decomposition (EEMD) technique to overcome the limitations of EMD by reducing the effects of noise and mode mixing in the decomposition process.However, EEMD has some drawbacks, such as increased computational complexity and the possibility of over-smoothing of the IMFs.Complete Ensemble Empirical Mode Decomposition with Adaptive Noise (CEEMDAN) technique 19 was introduced as an improvement to EEMD, where a fixed noise level is added to the signal to generate the ensemble realizations, CEEMDAN employs a noise level that incrementally grows and adjusts according to the specific characteristics of the local signal.CEEMDAN has another advantage over EEMD in that it is more effective at addressing the issue of mode mixing.Zhang and Yang 20 coupled CEEMDAN with Gated recurrent units (GRU) to forecast suspended sediment concentration in the Yellow River and outperformed Support Vector Machines, LSTM and GRU standalone models.Lu and Ma 21 suggested using CEEMDAN approach with XGBoost and Random Forest.Zhang et al. 22 suggested the decomposition of water quality data into IMFs using CEEMDAN algorithm and combining it with LSTM model for forecasting the water quality.Song and Yao 23 suggests using CEEMDAN decomposition in conjunction with LSTM for water quality prediction. 24ecommends using CEEMDAN and Variational Mode Decomposition (VMD) combined with Least Square Support Vector Machine (LSSVM) and Extreme Learning Machines (ELM) for estimating water quality parameters.
The ability of CEEMDAN to decompose the time series into IMFs facilitates the separation and analysis of various frequency components.This can help to detect the non-linear and non-stationary properties of the data.Deep learning and machine learning models have the capacity to learn and represent complex relationships between input variables, enabling them to capture both linear and non-linear dependencies in the data.This is especially useful in forecasting water quality since the interactions between input factors and water quality metrics can be complicated and multifaceted.This motivated us to investigate the potential of CEEMDAN decomposition in combination with AdaBoost and deep learning models in order to capture the long-term relationships accurately for DO forecasting.First, the CEEMDAN technique is used to decompose the data into a set of IMF sub-series components.The high and medium-frequency sub-series are forecasted using the Adaboost-BiLSTM, while the low-frequency components are forecasted using the standalone LSTM model.For a comprehensive analysis, the CEEMDAN-AdaBoost-BiLSTM-LSTM model is evaluated using data collected from ten different stations of river Ganga and compared with six standalone models and four models utilizing decomposition.Diebold Mariano test and two-sided t-test further suggests that there is a statistically significant difference in the forecasting performance of our proposed model as compared to the alternative models.

Study area and data collection
River Ganga, also known as the Ganges, is a transboundary river that flows through India and Bangladesh.It begins in the Himalayas, travels about 2525 kilometers, and empties into the Bay of Bengal.Apart from its cultural and spiritual significance, the Ganga River enables inland navigation and supports key sectors like as agriculture, fishing, tourism, and hydropower generation, all of which contribute to the local economy.Unfortunately, the water in the Ganga river is severely polluted, and the primary sources of which can be attributed to the unregulated release of untreated wastewater, industrial waste, and the inflow of agricultural runoff 25,26 .
River Ganga has a length of around 1450 kilometers in Uttar Pradesh.The data used in this study was collected from the Uttar Pradesh Pollution Control Board, Government of India.This dataset comprises of realtime measurements of DO levels obtained through the Real-time Water Quality Monitoring System installed at multiple locations along the River Ganga, its principal tributaries, and drains.The water quality data utilised in this study is for ten stations of River Ganga and its tributaries in Uttar Pradesh and spans from 1 April 2017 to 30 September 2021.In the original dataset, the sampling frequency ranged from one sample per hour for the years 2017 and 2018 to one sample every 15 min from the year 2019 to 2021.The specific locations analyzed in this study are distinctly marked and shown in Fig. 1.More details about the data considered in the study are provided in Supplementary Information and Supplementary Table S1.

CEEMDAN
CEEMDAN 19 uses an adaptive noise level, which adjusts the noise level of each ensemble realization based on the local signal characteristics, to enhance decomposition performance and efficiency.For the CEEMDAN decomposition process, first, we add white noise to the original time series signal X(t).This can be expressed as follows: where s 0 is the noise coefficient for controlling the signal-to-noise ratio and ω i (t) denotes the i th noise for i = 1...N , where N is the number of times EMD is performed.We use EMD to decompose each X i (t) into IMF i (t) .The first IMF component can be obtained by averaging all the modes.This can be mathematically expressed as follows: (1) www.nature.com/scientificreports/The initial residual is computed by subtracting the first IMF from the original signal as follows: The rest of the IMFs and residual can be calculated as follows: where, E k (.) extracts the k th IMF that EMD has decomposed.The calculation of the final residual involves sub- tracting the sum of all the K IMFs from the original signal, resulting in the following expression

AdaBoost
AdaBoost is a machine learning algorithm developed by by Freund et al. 27 .It employs a combination of weak classifiers to construct a strong classifier with enhanced performance.Although it was originally designed for classification problems, it can also be adapted for regression tasks 28 .In AdaBoost regression, weak learners are replaced with regression models that can predict continuous values.The algorithm iteratively fits regression models to the data and assigns weights to the samples based on their error.More weight is given to poorly predicted samples in subsequent iterations.The final prediction is a weighted combination of all the regression models' predictions, with the weights determined by their performance on the training data.AdaBoost uses the following steps for the computation: Let the data points used for training be (x 1 , y 1 ) ... (x N , y N ).
Step 1 Initialize the number of iterations as T. Let for the initial iteration, t = 1 and average loss L t = 0 .Initial- ize the sample weights distribution as D t (i) = 1 N for i = 1...N Step 2 Using the sample weights, train the regression model.Let f (x i ) be the predicted value for the i th sample.
Step 3 For each training sample compute the loss as : The loss function for each individual training instance can be calculated as: where, M t = max N i=1 l t (i) Step 5 Compute the average loss: www.nature.com/scientificreports/ Step 6 Set Step 7 Now, update the weight distribution of all the samples as follows : Here, Z t is the normalization factor.
Step 8 Now update t = t + 1 and repeat steps 2 to 8 while t ≤ T and average loss function L t ≤ 0.5 Step 9 Finally, the output can be given as follows:

LSTM
LSTM 29 are a special type of recurrent neural network with the capacity to learn long-term dependencies by employing cell states to store information about various time periods.A individual LSTM unit's cell state describes the data that has been thought to be relevant up to that date.By using input, output, and forget gates, LSTMs control the information flow to the cell states.The input gate assesses the importance of incoming data and decides how much of it should be retained in the memory cell.The forget gate regulates which past information should be ignored, while the output gate guarantees that only relevant information is produced at each time step.The input at the current time step, denoted as x t , and the previous hidden state, denoted as h t−1 , are used to calculate the values for the input gate, output gate, and forget gate.These gate values are then utilized to update the cell state c t of the LSTM.The output gate is applied to the hyperbolic tangent of the cell state c t to produce the final hidden state h t .The equations describing the working of LSTM are as follows: where, σ represents the sigmoid function, tanh denotes the hyperbolic tangent function and × indicates the element-wise multiplication.

BiLSTM
The BiLSTM layers are comprised of a pair of LSTM layers with one of them handling the input sequence in a forward direction and the other one processing it in reverse direction.BiLSTM can capture dependencies and patterns that a unidirectional LSTM would miss by analysing the input sequence in both directions resulting in a more complete contextual understanding.The equations describing the working of BiLSTM can be expressed mathematically as below: here, W fx and W bx are the weight matrices from the input to the recurrent units, and W fh and W bh are the weight matrices from the recurrent units to themselves for the forward and backward layers respectively.The biases for the forward and backward layers are given by b fb and b b .φ is the activation function at the hidden layers.W fy and W by are the weight matrices, and b y denotes the bias for the output layer.

Proposed approach
The steps for developing the model based on CEEMDAN decomposition combined with AdaBoost and deep learning are given below and the flowchart of the proposed model is shown in Fig. 2. where, x n corresponds to normalized value of x. x min and x max denotes the minimum and maximum value of the variable.5.For each IMF, Autocorrelation Function (ACF) and Partial Autocorrelation Function (PACF) are utilized to choose significant lagged data for constructing the forecasting models.Sliding window based multiple input multiple output approach is used to generate a collection of input-output pairs from time series data, with an input window providing historical information while the output window representing the target values.6.The zero-crossing rate is calculated to check the fluctuation frequency of the data in each IMF and then categorize them into high, medium and low frequency components.If zero-crossing rate of IMF ≥ 0.01 then the IMF is categorized as having high/medium frequency otherwise it has low frequency 30 .In order to forecast high and medium frequency IMFs, the AdaBoost-BiLSTM model is utilized due to its ability to predict complex non-linear patterns.BiLSTM can handle complex sequential data and detect long-term relationships between inputs and outputs, and AdaBoost can improve prediction accuracy by combining several BiLSTMs.LSTM models are employed for the low-frequency IMFs and residual that display reasonably simple and smooth patterns.In the end, the final prediction is obtained by summing the predicted outcomes of all the IMFs and residual.

Implementation details
The implemented model is developed using Python version 3.7.The CEEMDAN method is used with the help of EMD-signal 1.2.2 package.While implementing CEEMDAN method, the parameter settings include 100 trials and a scale of 0.05 for added noise.The number of lags needed as input to the forecasting algorithms is determined based on ACF and PACF.Supplementary Table S2 presents the identified significant time lags for the decomposed components of the time series data across ten stations.Each dataset may produce different number of IMFs on decomposition.Here, 'NA' as the table entry represents the absence of IMF for the particular station.Optimal hyperparameters are selected while building the models (details in Supplementary Information).The Supplementary Table S3 lists the search space explored for identifying the optimal hyperparameters.
(18) The zero-crossing rate of the IMFs for all the station datasets considered are given in Supplementary Table S4.
The deep learning models were implemented using TensorFlow 2.0.0 and Keras 2.3.1 packages.The number of hidden layers range from 3 to 5 and the number of neurons in the hidden layers is varied in the range of 8 to 64.
To improve the gradient descent technique, the Adam optimizer was utilized.The ReLU activation function is used in all the hidden layers of the models.The model is trained using the learning rate of 0.001.The batch size used for training is set to 64, and a variable number of epochs ranging from 10 to 100 are employed for each IMF.Early stopping is used to prevent overfitting of the model during training.

Performance evaluation metrics
The criteria to assess the model's performance are Mean Absolute Error (MAE), Mean Absolute Percentage Error (MAPE), Root Mean Square Error (RMSE) and Coefficient of Determination ( R 2 ) (details in Supplementary Information), mathematically, expressed as follows: where, n is the total count of the samples, y is the mean value, y t is the true value and ŷt is the predicted value at time t.

Results
Figure 3 shows the IMFs and residual that were generated using CEEMDAN approach for the DO time series data of Bithoor station.The time series is decomposed into 12 IMFs and a residual subseries.It can be seen that the frequency of first IMF is the highest, and as we progress toward the subsequent IMFs, the frequency gradually decreases.Based on zero-crossing rate, high and medium frequency IMFs are from IMF1 to IMF7, whereas IMF8 to IMF12 and residual are low frequency components.
We have used six single models without decomposition: Linear Regression (LR), SVR, Random Forest (RF), ANN, LSTM and BiLSTM.The models utilizing decomposition are: CEEMDAN-ANN, CEEMDAN-LSTM, CEEMDAN-BiLSTM, CEEMDAN-AdaBoost-BiLSTM and CEEMDAN-AdaBoost-BiLSTM-LSTM.The observed and forecasted DO values for all stations for the test data are presented in Fig. 4. The actual values are shown in black color, and blue, red, and green lines represent the predicted values for one, two and three-hour ahead forecast, respectively.Similarly, Fig. 5 illustrates the scatter plots showcasing the relationship between the observed and predicted DO values on test data.This corresponds to higher correlation between the measured and the predicted values.Tables 1 and 2 illustrates the comprehensive comparison of the individual models and hybrid models utilizing decomposition, respectively, highlighting their performance based on RMSE and MAE values.Figure 6 displays spider plots that compare the coefficient of determination values and help to visualize the significant improvement in the performance of CEEMDAN-AdaBoost-BiLSTM-LSTM over individual models Figure 3. IMFs and residual generated after performing CEEMDAN decomposition of the time series data at Bithoor station.for various stations.Figure 7 illustrates the MAPE values comparing all methods with the proposed approach for all horizons for every station.It can be seen that the MAPE value for the proposed approach is lower than all other models.This demonstrates that the proposed approach has the potential to serve as an efficient means of generating dependable and accurate forecasting.

Discussion
Among the individual models, it can be observed that the LR model gives satisfactory performance in several cases.However, it is important to note that, LR is more capable to capture linear relationships within the data.On the other hand, SVR and RF are capable of handling the non-linear and complex relationships.The ANN and deep learning models like LSTM and BiLSTM are proficient in capturing non-linear and long-term temporal dependencies in DO time series data.In terms of all evaluation metrics, the BiLSTM model consistently demonstrates superior forecasting accuracy across all forecasting horizons.
By utilizing the CEEMDAN technique to decompose the water quality data into various frequency segments and subsequently applying the deep learning models to capture the temporal dependencies within each frequency segments, a substantial improvement in forecasting accuracy has been observed as compared to the models only utilizing standalone deep learning model.The effectiveness of the CEEMDAN decomposition   The CEEMDAN-AdaBoost-BiLSTM model demonstrates up to 15.517%, 12.290% and 11.111% reductions over CEEMDAN-BiLSTM for one, two, and three-hour ahead predictions respectively for RMSE, across all stations.Similarly, for MAE reductions of up to 21.519%, 13.803%, and 18.182% across all stations is noticed for one, two, and three-hour ahead predictions respectively.This suggests that using AdaBoost can enhance the performance when compared with using only an individual BiLSTM model.Training multiple BiLSTM models in conjunction with the Adaboost algorithm involves a dynamic adjustment of model weights based on prediction errors.This collaborative approach serves to enhance forecasting outcome and robustness.
However, using the same approach to model all the components might not be optimal.The different IMF components are likely to exhibit distinct characteristics.High and medium-frequency IMFs, marked by swift fluctuations and intricate patterns, can be accurately predicted with AdaBoost-BiLSTM.These IMFs signify short-term variations and exhibit high randomness.Conversely, low-frequency IMFs and residual components, 323% are seen across one, two, and three-hour ahead predictions respectively.This approach acknowledges the variability among components and emphasizes the importance of adapting prediction techniques to the individual characteristics of each component in order to produce more accurate results.The proposed approach has consistently provided the best forecast across all forecast horizons and all stations, thus demonstrating its superiority.It can be observed that, as the prediction time step increases, the predictive performance of all models gradually declines, demonstrating an increasing accumulation of errors in multi-step forecasting.This phenomena correctly depicts that predicting further into the future gets more difficult as the time step lengthens.However, improved outcomes are still obtained for our proposed CEEMDAN-AdaBoost-BiLSTM-LSTM model.Although decomposing the original data is an additional computational step in the modeling process, it significantly enhances prediction accuracy.The simple and less complex decomposition-based models apply the same technique to forecast all the IMFs without considering their characteristics.While this approach may seem straightforward, it comes at the expense of predictive accuracy.The proposed CEEMDAN-AdaBoost-BiLSTM-LSTM model introduces an additional step of computing the zero-crossing rate for each IMF, aiding in selecting the appropriate forecasting model (AdaBoost-BiLSTM or LSTM) for each IMF, thereby improving accuracy.This computation incurs only a minimal increase in computation overhead and execution time, taking approximately 0.0114 seconds (example of Bithoor dataset).The application of AdaBoost-BiLSTM to the complex high and medium frequency IMFs further enhances predictive performance through ensemble learning techniques, making it a useful approach in critical real-world applications like water quality monitoring and forecasting.
Diebold-Mariano (DM) test 31 and t-test are commonly used statistical tests used to compare the performance of models relative to each other [32][33][34][35] .DM test with the squared error as the loss function is used to compare the proposed forecasting model with ten alternative models.The DM statistics obtained are compiled in Supplementary Table S5.First, the magnitude of the DM value shows that the proposed model has a significant advantage over the other models in terms of forecasting accuracy.Furthermore, the computed p-values play a critical role in determining the statistical significance of the observed differences.Notably, each p-value was found below the 0.05 significance level.This suggests that the forecast accuracy of the proposed model differs statistically by a significant amount from that of the alternative models.Additionally, we conducted two-sided t-tests using results from seven different executions using different random seeding to evaluate the significance of performance variations among different datasets, as outlined in Supplementary Table S6.Notably, the p-values consistently remained below the 0.05 significance level, indicating that the proposed model's RMSE scores significantly outperformed those of the benchmark models across all datasets.
Considering the dataset for the Bithoor station, the average runtimes for the individual models are as follows: 7.35 seconds for LR, 15.33 seconds for SVR, 10.32 seconds for RF, 90.20 seconds for ANN, 117.23 seconds for LSTM and 126.11 seconds for BiLSTM.The models with decomposition as a pre-processing step have the following average runtimes: 910.02 seconds for CEEMDAN-ANN, 1233.71seconds for CEEMDAN-LSTM, 1254.75 seconds for CEEMDAN-BiLSTM, 3143.23 seconds for CEEMDAN-AdaBoost-BiLSTM and 1656.70 seconds for CEEMDAN-AdaBoost-BiLSTM-LSTM.The data decomposition pre-processing step increases the runtime of the hybrid models.Out of all the models, CEEMDAN-AdaBoost-BiLSTM takes the longest to run.It is worth noting that while the CEEMDAN-AdaBoost-BiLSTM model may show comparable accuracy levels to our proposed CEEMDAN-AdaBoost-BiLSTM-LSTM model in a few isolated cases, it does so at the cost of increased execution time.This is because the former model applies AdaBoost-BiLSTM to all decomposed components, whereas our proposed model applies AdaBoost-BiLSTM only to the high and medium frequency IMFs.By applying AdaBoost only to specific components, our model significantly reduces the execution time.

Conclusion
This study examines the use of a novel approach combining the CEEMDAN approach in conjunction with AdaBoost and deep-learning for short-term and multi-step forecasting of the DO levels in river Ganga.The experiments demonstrate that: 1.A more precise and superior prediction result is achieved by taking into account the unique qualities of each IMF component and using a variety of prediction techniques catered to these components.Here, the AdaBoost-BiLSTM model is used to predict the high and medium frequency IMFs with complex non-linear patterns and standalone LSTM model is used for the low-frequency IMFs and residual that display reasonably simple and smooth patterns.The proposed approach outperforms all the models that used the same technique to forecast each component.

Figure 1 .
Figure 1.The geographic locations of the ten stations of River Ganga in Uttar Pradesh, India are denoted by blue markers.
W f , W i , W o , W c are the input weight matrices and U f , U i , U o , U c are the recurrent weight matrices for the forget, input, output and memory cell gate respectively.The bias vectors to the respective gates are represented by b f , b i , b o and b c .The hidden state and input at timestamp t are denoted by h t and x t respectively.
) ŷt =σ (W fy h ft + W by h bt + b y ) Vol.:(0123456789) Scientific Reports | (2024) 14:11199 | https://doi.org/10.1038/s41598-024-61910-w 1.The first step is to use linear interpolation to fill in the missing values.The data for all the years are combined.To ensure consistent sampling frequency across the dataset, the data was extracted at a constant sample rate of one observation per hour for all years from the combined dataset.2. The CEEMDAN approach is used to decompose the DO time series data into several IMFs and residual with each generated IMF possessing unique and distinctive inherent properties.3. The data is partitioned into training and testing sets with a ratio of 75:25.10% of the training data is taken for validation.4. The data is normalized using min-max normalization, which scales the values to the [0,1] range as follows:

Figure 2 .
Figure 2. Flowchart of the proposed approach.

Figure 4 .
Figure 4.A visual representation of actual and predicted DO values over the test data at different stations.

Figure 5 .
Figure 5. Scatter plot illustrating the alignment between observed and predicted DO values on test data.

Figure 7 .
Figure 7.Comparison of MAPE values for different models over test data across all stations.

2 .
The proposed CEEMDAN-AdaBoost-BiLSTM-LSTM model outperforms the CEEMDAN-AdaBoost-BiL-STM with reductions in RMSE by up to 27.491%, 23.280%, and 11.567% and the reductions in MAE by up to 31.122%, 34.317%, and 21.323% for 1, 2, and 3-hour ahead predictions respectively.Compared with CEEMDAN-BiLSTM, the proposed model demonstrates up to 27.986%, 30.622% and 17.767% reductions for RMSE, and up to 38.356%, 39.661%, and 34.154% reductions in MAE for 1, 2, and 3-hour ahead predictions respectively.Results of Diebold-Mariano and t-test suggests statistically significant difference in forecasts between our proposed and the models used in comparison.3. The proposed CEEMDAN-AdaBoost-BiLSTM-LSTM model is comparatively more computationally efficient than CEEMDAN-AdaBoost-BiLSTM model.The proposed model applies AdaBoost-BiLSTM only to the high frequency and medium frequency IMFs and employs a standalone LSTM model to the low frequency IMFs and residual.This takes lesser time to execute than CEEMDAN-AdaBoost-BiLSTM which applies AdaBoost-BiLSTM to all the IMFs and residual.However, the proposed model requires slightly longer processing time compared to the other simpler models.The future research will focus on integrating parallelization techniques for decreasing the computational time of the model making it more practically applicable while optimizing the performance.
Likewise, for CEEMDAN-LSTM model, the reductions in RMSE range from 25.556% to 54.229%, 13.610% to 50%, and 7.043% to 49.315%, for one-hour, two-hour, and three-hour ahead forecasts respectively, compared to the simple LSTM model.In terms of MAE, the reductions range from 5.063% to 45.544% for the one-hour ahead forecasts, 7.143% to 58.871% for the two-hour ahead forecasts, and 4.640% to 49.514% for the three-hour ahead forecasts.The CEEMDAN-BiLSTM exhibits superior performance by giving lower RMSE values compared to the simple BiLSTM, with reductions ranging from 22.661% to 54.271%, 7.335% to 52.877% and 3.895% to 48.732% for one, two, and three-hour ahead predictions across stations.Likewise, for MAE, reductions in the range of 1.357% to 47.535%, 3.488% to 49.484%, and 3.325% to 48.979% are observed for one-hour, two-hour, and three-hour forecasts, respectively.These results indicate that the approaches using decomposition as a data preprocessing step has a substantial influence on the efficacy of models.
org/10.1038/s41598-024-61910-w www.nature.com/scientificreports/method in improving the precision of water quality forecasts is evident on comparing the RMSE and MAE values.When compared to standalone ANN, the CEEMDAN-ANN demonstrates better performance by achieving lower RMSE values.The reductions vary from 8.849% to 55.656%, 10.954% to 51.213%, and 7.374% to 48.359% across all stations for one, two, and three-hour ahead forecasting horizons, respectively.In terms of MAE, the reductions vary from 2.109% to 60.811%, 0.414% to 52.290%, and 10.782% to 51.440% for one-hour, two-hour, and three-hour ahead forecasts.

Table 1 .
Prediction results on test data for single models without decomposition.

Table 2 .
Prediction results on test data for decomposition based hybrid models.Comparison of coefficient of determination values for different models over test data across all stations.primarilyrepresenting gradual trends or periodic elements, are more effectively predicted using standalone LSTM.The proposed CEEMDAN-AdaBoost-BiLSTM-LSTM model outperforms CEEMDAN-AdaBoost-BiL-STM with reductions of up to 27.491%, 20.375%, and 11.567% for one, two, and three-hour ahead predictions respectively for RMSE, across all stations.For MAE, the reductions up to 31.122%, 34.317%, and 21.